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Abstract. We presented the results of several statistical tests of the randomness in the angular 
sky-distribution of gamma-ray bursts in BATSE Catalog. Thirteen different tests were presented 
based on Voronoi tesselation, Minimal spanning tree and Multifractal spectrum for five classes 
(short 1, short2, intermediate, longl, long2) of gamma-ray bursts, separately. The longl and long2 
classes are distributed randomly. The intermediate subclass, in accordance with the earlier results of 
the authors, is distributed non-randomly. Concerning the short subclass earlier statistical tests also 
suggested some departure from the random distribution, but not on a high enough confidence level. 
The new tests presented in this article suggest also non-randomness here. 
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INTRODUCTION 

There are increasing evidences that all the GRBs do not represent a physically homo- 
geneous group (Kouveliotou et al. 1993, horvath 1998, balazs et al. 2003, Hakkila et al. 
2003, horvath et al. 2006). Hence, it is worth investigating that the physically different 
subgroups are also different in their angular distributions. In the last years the authors 
provided (balazs et al. 1998, balazs et al. 1999, Meszaros et al. 2000) several different 
tests probing the intrinsic isotropy in the angular sky-distribution of GRBs collected in 
BATSE Catalogs. One may conclude the results of these studies: A. The long subgroup 
seems to be distributed isotropically (see also Briggs 1993); B. The intermediate sub- 
group (horvath 2002, horvath et al. 2006) is distributed anisotropically on the 96-97% 
significance level; C. For the short subgroup the assumption of isotropy is rejected only 
on the 92% significance level; D. The long and the short subclasses, respectively, are dis- 
tributed differently on the 99.3% significance level. Independently (Litvin et al. 2001), 
confirmed the results A., B. and C. with one essential difference: for the intermediate 
subclass a much higher - namely 99.89% - significance level of anisotropy is claimed. 
Again, the short subgroup is found to be "suspicious", but only on the 85-95% signifi- 
cance level. In this paper, similarly to the previous studies, the intrinsic randomness is 
tested; this means that the non-uniform sky-exposure function of BATSE instrument was 
considered. 



TABLE 1. Tested samples of BATSE GRBs. 



Sample 


Duration 


Peak flux 


Number 




[s] 


[photons cm~^s~^] 


of GRBs 


Shortl 


Tgo < 2 s 


0.65 < ^256 < 2 


261 


Short2 


JgO < 2 S 


0.65 < ^256 


406 


Intermediate 


2 S < Tgo < 10 s 


0.65 < P256 


253 


Longl 


Tgo > 2 s 


0.65 < ^256 < 2 


676 


Long2 


T90 > 10 s 


0.65 < ^256 


966 



MATHEMATICAL SUMMARY AND THE TEST- VARIABLES 

The randomness of the point field on the sphere can be tested with respect to different 
criteria. In the following we defined several test-variables. 

Voronoi tesselation (VT). The Voronoi diagram - also known as Dirichlet tesselation 
or Thiessen polygons - is a fundamental structure in computational geometry (Voronoi 
1908, Stoyan & Stoyan 1994). Generally, this diagram provides a partition of a point 
pattern according to its spatial structure. Assume that there are N points (N » 1) scattered 
on a sphere surface with an unit radius. The Voronoi cell of a point is the region of the 
sphere surface consisting of points which are closer to this given point than to any other 
ones of the sphere. This cell forms a polygon on this sphere. Every such cell has its 
area (A) given in steradians, perimeter (P) given by the length of boundary (one great 
circle of the boundary curve is called also as "chord"), number of vertices (A'v) given 
by an integer positive number, and by the inner angles. This method is completely non- 
parametric, and therefore may be sensitive for various point pattern structures in the 
different subclasses of GRBs. 

Any of the four quantities characterizing the Voronoi cell can be used as test- variables 
or even some of their combinations, too. We defined the following test- variables: 1, 
Cell area A; 2, Cell vertex (edge) A'v; 3, Cell chords C; 4, Inner angle a,; 5, Round 
factor (RF) average RFav = AnA/P\ 6, Round factor (RF) homogeneity 1 — 7, 

Shape factor A/P^; 8, Modal factor o{ai)/Nv\ 9, The so-called "AD factor" defined as 
AD=l-(l-o{A)/{A))-\ 

Minimal spanning tree (MST). Contrary to VT, this method considers the distances 
(edges) among the points (vertices). A spanning tree is a system of lines connecting 
all the points without any loops. The minimal spanning tree (MST) is a system of 
connecting lines, where the sum of the lengths is minimal among all the possible 
connections between the points (Prim 1957). The statistics of the lengths and the MST 
angles between the edges at the vertices can be used for testing the randomness of the 
point pattern. To characterize the stochastic properties of a point patters we use three 
quantities obtained from a MST: 1, Variance of the MST edge-length o{Lmst) ; 2, Mean 
MST edge-length Lmst\ 3, Mean angle between edges olmst- 

Multifractal spectrum is the third method which was used. Here the only used vari- 
able is the f{a) multifractal spectrum, which is a sensitive tool for testing the non- 
randomness of a point pattern. 



TABLE 2. Calculated significance levels for the 13 test- variables and the five samples. 
A significance greater than 95% is put in bold face. 



Name 


var 


shortl 


short2 


interm. 


longl 


long2 


Cell area 


A 


36.82 


29.85 


94.53 


79.60 


82.59 


Cell vertex (edge) 


Nv 


36.82 


87.06 


2.99 


26.87 


7.96 


Cell chords 


C 


47.26 


52.24 


18.91 


84.58 


54.23 


Inner angle 


«, 


96.52 


21.39 


87.56 


37.81 


63.18 


RF average 


A%A/P 


65.17 


99.98 


33.83 


10.95 


86.07 


RF homogeneity 


^ RFav 


19.90 


A 'JO 

24.38 


CO T 1 

58.71 


55.72 


OA 

32.84 


Shape factor 




91.04 


94.03 


90.05 


55.22 


63.68 


Modal factor 


a{a,)lN, 


97.51 


1.99 


7.46 


56.22 


8.96 


AD factor 




32.84 


25.37 


11.44 


95.52 


52.74 


MST variance 




52.74 


38.31 


22.39 


13.93 


59.70 


MST average 


Lmst 


97.51 


7.46 


89.05 


56.72 


8.96 


MST angle 




85.07 


14.43 


36.82 


73.63 


60.70 


MFR spec Ira 


.f(cc) 


95.52 


96.02 


98.01 


73.63 


36.32 


Binomial test 




99.79 


99.74 


77.00 


55.13 




Squared EucUdean 


distance 


99.90 


99.98 


98.51 


93.03 


36.81 



RESULTS 

Completing 200 simulations in all of the subsamples (for them see Table 1.) we get a 
13D sample representing the joint probability distribution of the 13 test variables. Using 
a suitable chosen measure of distance of the points from the sample mean we can get 
a stochastic variable characterizing the deviation of the simulated points from the mean 
only by chance. An obvious choice would be the squared Euclidean distance. 

In case of a Gaussian distribution with unit variances and without correlations this 
would resulted in a distribution of 13 degree of freedom. But the test variables in our 
case are correlated and have different scales. Factor analysis (FA) is a suitable way to 
represent the correlated observed variables with fewer non-correlated variables of less 
in number (Wallet & Dussert 1998). The number of non-correlated variables, k, can be 
constrained by k < 8.377 in our case for n = 13. Hence, we retained 8 non-correlated 
variables. 

Out of the 13 test-variables only the multifractal spectrum gave significant (>95%) de- 
viation from the simulated sample in more than one group. The BATSE samples, how- 
ever, were different in the number of test-variables giving positive signal (>95%) and 
in the level of significance. Among the tested samples shortl experiences four (96.5%, 
97.5%, 97.5%, 95.5%), short2 two (99.98%, 96.02%), intermediate one (98.0%), longl 
one (95.5%) and long2 no variables with >95% significance (see Table 2.). Calculating 
the joint significance level we assumed that they can be represented as a linear combina- 
tion of non-correlated hidden factors of less in number. We obtained k=8 as the number 
of hidden factors. Then we computed the distribution of the squared Euclidean distances 
from the mean of the simulated variables. Comparing the distribution of the squared dis- 
tances of the simulated with the BATSE samples we concluded that the shortl, short2 
and intermediate groups deviate significantly (99.90%, 99.98% and 98.51%) from the 
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FIGURE 1. Distribution of the Euclidean distances of the simulated samples from the stochastic mean 
of the variables in the 13D parameter space. There are altogether 1000 simulated points. Full line marks 
a distribution of 8 degree of freedom, normalized to the sample size. The distances of the BATSE 
samples are also indicated. The departures of samples "shorti" and "short!" exceed all those of the 
simulated points. The probabilities, that these deviations are non-random, equal 99.9% and 99.98%. 



fully randomness but it is not the case at the long samples (see Fig. 1.). 

Acknowledgments: This study was supported by OTKA grant No. T048870 and 75072 
and Res. Program MSM0021620860 and GAUK grant No.46307 (A.M.). 



1. Balazs E.G., Meszaros A., Horvath I., 1998, A&A, 339, 1 

2. Balazs E.G., Meszaros A., Horvath 1., Vavrek R., 1999, A&A Suppl., 138, 417 

3. Balazs E.G., Bagoly Z., Horvath 1., Meszaros A., Meszaros P., 2003, A&A, 401, 129 

4. Briggs M., 1993, ApJ, 407, 126 

5. Hakkila, J., et al., 2003, ApJ, 582, 320 

6. Horvath E, 1998, ApJ, 508, 757 

7. Horvath 1., 2002, A&A, 392, 791 

8. Horvath 1., Balazs E.G., Bagoly Z., Ryde P., Meszaros A., 2006, A&A, 447, 23 

9. Kouveliotou, C., et al., 1993, ApJ, 413, ElOl 

10. Eitvin V.F., Matveev S.A., Mamedov S.V., Orlov V.V., 2001, Pis'ma v Astron. Zhurnal, 27, 489 

1 1 . Meszaros A., Bagoly Z., Horvath E, Balazs E.G., Vavrek R., 2000b, ApJ, 539, 98 

12. Prim R.C., 1957, Bell Syst. Techn. Journ., 36, 1389 

13. Stoyan D., Stoyan H., 1994, Fractals, Random Shapes and Point Fields, Wiley J. & Sons, New York 

14. Voronoi G., 1908, J. Reine Angew. Math., 134, 198 

15. Wallet, F., Dussert C., 1998, Europhys. Eet., 42, 493 



REFERENCES 



